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, A description of the transition from the inflationary epoch to radiation 

domination requires the understanding of quantum fields out of thermal 
equilibrium, particle creation and thermalisation. This can be studied 
from first principles by solving a set of truncated real-time Schwinger- 
5^ \ Dyson equations, written in terms of the mean field (inflaton) and the field 

propagators, derived from the two-particle irreducible effective action. We 
investigate some aspects of this problem by considering the dynamics of 
a slow- rolling mean field coupled to a second quantum field, using a (p'^x^ 
interaction. We focus on thermal effects. It is found that interactions lead 
to an earlier end of slow-roll and that the evolution afterwards depends on 
details of the heatbath. 
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1 Introduction 



Cosmological observations strongly suggest that a stage of accelerated expansion 
took place in the very early Universe [U [2]. Such a stage is well described in 
terms of the dynamics of a scalar field, the inflaton = {(p(t,x.)), slow-rolling 
in a suitable potential V[(f)]. Classically, the inflaton equation of motion reads 

0(t)+3i7(t)0(t) + V'[0(t)] = O, (1.1) 

with the Hubble rate determined by the Friedmann equation 

i4,^t) + vm], 

^ = ^ 3MJ -' <1-^' 

where Mp\ is the Planck mass. 

It is natural to assume that the inflaton is coupled directly or indirectly to the 
Standard Model fields and that after inflation, the energy stored in the inflaton 
is transfered to excitations of these fields. This can be realised through pertur- 
bative reheating [3l H] or nonperturbative preheating O El [TJ [H [9l [TOl [HI [121 [IS] 
processes. Through scattering, the system eventually reaches thermal equilib- 
rium and a radiation dominated Universe emerges. In quantum field theory, the 
inflaton evolution in real time is described by a Schwinger-Dyson equation, 

0(t) + 3H{t)^{t) + V'[<j){t)] = - [ dt' S^(t, t')<f){t'), (1.3) 

J to 

where to = is the initial time, taken as the starting point for the evolution. The 
nonlocal (self-energy like) term Il^{t,t') contains the interaction with quantum 
fields and is to be determined, either through some well-motivated ansatz or in 
terms of a systematic diagrammatic expansion. 

Due to the interactions, one may expect a backreaction on the inflaton re- 
sulting in, e.g., a modified effective mass parameter and possibly friction. A 
phenomenological inflaton equation of motion incorporating this would be 

0(t) + [3Hit) + T(t)] 0(t) + Kff[0(t)] = 0. (1.4) 

Various approximations have been invoked to motivate Eq. fll.41) [TH [TJl [ISl |T71 
[18] F1 Effective inflaton dynamics as described by Eq. (11.41) is important for warm 
inflation [T5l [22] . in which it is assumed that particles created during inflation 
interact fast enough such that a quasi-stable nonvacuum state is reached. 

For a mean field 0(t) undergoing small amplitude oscillations in a thermal 
background, an effective damping rate T is indeed generated, and it can be cal- 
culated perturbatively in a linear response treatment [23]. However, when the 



^See also Refs. [T9l[20l[2T] for objections to this equation. 
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mean field displacement from equilibrium cannot be treated as a small pertur- 
bation, as in the case of a slow-rolling infiaton, things are less clear and a fully 
nonequilibrium treatment is required. This is even more the case if properties of 
the heatbath depend on the value of the mean field through interactions. 

In this paper we will not attempt a full investigation in an expanding back- 
ground. Instead, we will study some aspects of this problem in a simple inflation- 
inspired model using the tools of nonequilibrium field theory, i.e. solving the 
Schwinger-Dyson equation (11.31) . as derived from truncations of the two-particle 
irreducible (2PI) effective action [211 [25]. section [2] we introduce the model, ex- 
plain how we deal with the expansion and describe the dynamics of the mean field 
and fluctuations in absence of most interactions. We write down the Schwinger- 
Dyson equations to next-to-leading order in a coupling expansion of the 2PI effec- 
tive action in Section [31 In Section [4] we subsequently solve the lattice discretised 
integro-differential equations numerically without further approximations. We 
consider various scenarios and study in particular the role of thermal initial con- 
ditions for the two quantum fields in our model. Our findings are summarized in 
Section[5l In the Appendix we collect several approximations that are well-known 
in the literature and that are used in the main body of the paper for comparison. 



2 Model and free evolution 

We consider two interacting scalar fields if and x evolving in a flat Friedmann- 
Robertson- Walker Universe, with the action 



S = — / d X 



(2.1) 

The metric is specified by ds^ = dt^ — a?{t)dy? . Hence, the spatial derivatives 
are given by 9x/a(t), where x is co-moving. The Friedmann equation determines 
the evolution of the scale factor as 

-2 

pi 



, , , _ HHt) = ^ Y^, (2.2) 

a(t) J ^ ' 3M2 ' ^ ^ 



where (T°°) is the renormalized expectation value of the energy density of matter. 
We let if play the role of the infiaton: its expectation value {ip{t,x.)) = (pit) is 
homogeneous and set to some large value initially, 0(0) = 00, prompting slow-roll 
infiation. We take (x(^,x)) = for all times, consistent with the symmetries of 
the action. 

In order to understand the dynamics in this model, we start by considering 
the classical evolution of 0(t) without interactions. Its equation of motion is 

0(t) + 3^r(t)0(t)+mXt) = O, (2.3) 
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with H given by 



H\t) = ^ ^ ^ r ■ (2.4) 
In the slow-roll limit, H'^if) = m^0^(t)/6Mpi, and we find 

<Pit) = 00 - )J^m^Mp,t. (2.5) 

Assuming that infiation ends when the slow-roll parameter e = |Mp[(V^'[0]/V^[</)])^ 
equals 1, we find that there is infiation when 0/Mpi > \/2, irrespective of the 
value of m^; a given initial value 0o corresponds in the slow- roll approximaton 
to Ne = [{(j)o/Mpiy - 2]/4 e-folds before the end of infiation. 

The X field, coupled to this slow-rolling infiaton, is expanded in terms of 
creation and annihilation operators as 

X(t, x) = ^ {aMt)e'^'' + 4/k We"^*^-^) , (2.6) 



where 



(2.7) 



k J (2vr)3- 
The mode functions fkit) obey the equation 

[d^ + 3Hm + kV«'(^) + g'<P\t)] /k(t) = 0. (2.8) 

At this level the infiaton acts as a time-dependent mass. 

We now make a practical choice that we will follow throughout this paper: we 
ignore the expansion of the Universe for the quantum fields (or mode functions) 
and put H = 0, a{t) = 1 in Eq. (12.81) . while keeping Hubble friction in the 
mean field equation. This has the effect of emphasizing the role of scattering and 
energy transfer during the slow-roll regime and maximises the backreaction of the 
quantum modes on the mean field by omitting the effect of dilution and redshift. 
As a result we are not studying infiation, but rather the impact of interactions on 
a slow-rolling mean field. It allows us to disentangle the direct effect of field theory 
interactions from the effect of including a x-component in T^^. Such a component 
would influence the inflaton evolution via the Friedmann equation. Neglecting 
expansion is expected to strongly favour thermalisation. We can make a rough 
estimate by comparing the expansion rate H to the perturbative damping rate 
r(T) [23]. The ratio T/H is maximal at the end of inflation, and for the largest 
temperatures used in Section H] (T/Mpi = 4), we flnd F/if ~ 0.03. The expansion 
could therefore have a signiflcant impact on thermalisation. In practice however, 
the time scales under consideration here are much shorter than the thermalisation 
time scale, and, as we will see below, equilibration and thermalization will not 
play a role here. 
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Figure 1: The equal-time correlator {x^it, x)) — (x^(0, x)) for free evolution of the 
X modes in the mean-field background. Shown are the numerically determined 
evolution in the exact (full) and the slow-roll (dashed) background, as well as 
in the adiabatic approximation (long-dashed). The inset shows the mean field 
during the exact evolution (full) and in the slow-roll approximation (dashed). 

In Fig. [T] we show the evolution of the mean field (j){t) and the (subtracted) 
equal-time correlator 

(X^(t,x))= / \Mt)\\ (2.9) 

for an initial state in vacuum. The initial value of the mean field is (po/M^i = 
10 > V2, ensuring an extended slow-roll stage. We show the evolution of (x^) 
determined by numerically solving the equations for the mode functions fkit) 
using various approximations detailed in Appendices lA.ll and IA.2I in the exact 
mean-field background determined by Eqs. fl2.3|2.4p . in the slow-roll background 
(12. 5p . and in the adiabatic approximation. In the inset the time evolution of the 
mean field is shown, determined by solving Eqs. (12. 3112. 4p numerically and in the 
slow-roll approximation (12. 5p . In the latter the evolution is stopped when 0(t) 
reaches zero; (j){t) is kept constant at zero from then on. 

During the slow-roll regime the mean field (j){t) and hence the effective x 
mode energies ujk{t) = [k^ -|-m^ -|- (yf^0^(t)]^/^ decrease. As a result the equal-time 
correlator (x^) increases in time. The sharp increase is described very well in the 
adiabatic approximation. Therefore it does not correspond to a large amount of 
particle production. The question we address in the remainder of the paper is 
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whether and how interactions modify this basic scenario. 



3 2PI expansion at weak coupling 

Corrections to the dynamics considered above arise from quantum and possibly 
thermal back-reaction effects. These can be deduced from taking the expectation 
value of the Heisenberg operator equation of motion, yielding 

dl{^{x)) + 3H{t)dt{^{x)) + mli^ix)) + g'ix'ixMx)) = 0. (3.1) 

The last term can be expanded in terms of connected correlators, 

{X\XMX)) = {X\x))<j){t) + (x'(x)(/^(x))connected, (3.2) 

using again that {x{x)) = and writing {(p{x)) = (pit). As is well-known, keeping 
only the first term amounts to the Hartree or Gaussian approximation (see Ap- 
pendix [A]3]) . In that case the effects of interactions reduce to a time-dependent 
mass. In order to systematically go beyond the Hartree approximation and in- 
clude effects from the connected three-point correlator in Eq. (13.21) . we use a 
formalism that allows for consistent truncations of the Schwinger- Dyson hierar- 
chy: the two-particle irreducible (2PI) effective action. 

The 2PI effective action is a functional written in terms of the mean field 
and the propagators G^p, G^ 



r[0, G^, = S[4>] + ^Tr In G'i + ^Tr G^i^^.o 

+ ^IVlnG-i +^TrG-iG^,o + r2[0,G^,GJ, (3.3) 

where S'[0] is the tree-level action (12.11) written in terms of the mean field, G~j^ q 
are the free inverse propagators, 

zG^^^ix, y) = - + ml] 5c{x - y), (3.4) 
^Go,i(x, y) = - [D, + ml + /^'(x)] 5c{x - y), (3.5) 

and r2 is the sum of all two-particle-irreducible diagrams, depending on full 
propagators and the mean field via the three-point vertex g'^ipipy^. We follow 
the notation of Ref . [26] . 

The series of diagrams in r2 can be truncated, with each set of diagrams cor- 
responding to an approximation to the full theory. In a weak-coupling expansion 
we include (see Fig. |2]), 

2r2[0, G^, = 1" d'^xG^{x,x)G^{x,x) (3.6) 
4 Jc 

j d^xd'^y(l){x)G^{x,y)Gl{x,y)(l){y) -^-^ J d'^xd'^y Gl{x,y)Gl{x,y). 
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a) b) c) 



Figure 2: Diagrams included in the 2PI effective action: a) figure 8 b) pippi 
c) basketball. Full/dashed lines are xlv propagators. The crosses denote the 
inflaton mean field. 



Diagram b) is formally ~ (7^0^ • We have ignored the three-loop diagram that is 
~ g^c^y" . All g^c^'^ insertions are resummed already at the Hartree level: therefore 
the truncation we use here formally amounts to assuming that g^c^ is small, 
but with no constraints on (7^0^. Convergence properties of the 2PI coupling 
expansion were briefly discussed in Ref. [23], and for the expansion in the 
0{N) model in Refs. [271 [28]. 

The equations of motion result from taking the functional derivatives. 



5r [0,G'^,GJ 



5r[0,G^,G,] ^^^^^ 



6G^{x,y) 



After going to momentum space and writing the contour propagators and self 
energies in terms of the statistical and spectral components. 



G^{t, t'- k) = t'- k) - - signc(t, t')p^{t, t'- k). 



we find the standard equations [251 ESI ED] 

[d^, + k^ + Mli^it)] p^/^it, t'- k) = - [ dt" S^/^,p(t, t"- k)p^/^(t", t'- k), 

Jt' 

+ k^ + Mj/^(t)] F^/^it, t'; k) = - [ dt" S^/^,,(t, t"; k)F^/^{t", t'; k) 

Jo 

+ [ dt"E^/^,Fit,t";k)p^/^{t",t';k), 



(3.9) 



with 



Ml{t)=ml + g' jj^{t,t-k), 



M^{t)=ml + g' 



(j)'{t)+ / F^{t,t;k) 



(3.10) 
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The nonlocal self-energies originate from diagrams b) and c) and are easiest 
written in real space. They are 



(3.11) 



where 



^xMF*^^'^) = -2/0(a;) [F'xi^^y)K/xi^^y) - p'fi/xi^^y)pxi^^y)/^] ^iv)^ 



and 



- V^^/x(a;, i/)p</./x(a;' i/)^x/^(3;' y) 



K/xA""^ y) = -2^?' 2/) - Px/^(^' ^)/4] F^/xi^, y) 

+9^Fx/A^^ y)px/'pi^^ y)pv/xi^^ y)^ 



The mean field equation is given by 



with Ml{t) = M^{t), H{t) determined by Eq. ([M]), and 

E^(t,t') = Sa(t,t';k = 0). 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



4 Numerical analysis 

In order to solve the set of coupled integro-differential equations presented in 
the previous sections, the system is discretized on a space-time lattice with spa- 
tial lattice spacing and temporal spacing at- The resulting discretized 2PI 
equations are solved numerically, see Refs. [30l [25| [23] and references therein for 
details. The initial density matrix is taken to be Gaussian, such that only the 
one- and two-point functions need to be initialized. For the one-point functions 
we take 0(0) = but a large initial 0(0) = 00 = lOMpi. For the two-point 
functions we use 



G^/^(t,t'; k) 
dtdt'G^/^{t,t'; k) 



n 



x/vfi 



+ 1/2 



t=t'=0 



t=t'=0 



t=t'=0 



UJ- 



x/v^,o 



n 



x/'Pfi 



+ 1/2 cu, 



,X/'P,0 



(4.1) 

(4.2) 
(4.3) 
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where the initial mode energy is 



k2 + M2 (0), (4.4) 



and the initial particle number is 



= [ exp 



^k'''"'7T,/, -1) . (4.5) 



^1 



For the lattice discretized formulation, an 0{at) improvement term is added to 
the RHS of Eq. (gj]), 



1 [G,/,(a„ 0; k) - G,/,(0, 0; k)] = -| [n^J-'' + I/2) u^J^'', (4.6) 



which ensures that, in the Hartree approximation, evolution initialized at the 
Hartree fixed point [21] remains at that fixed point, for finite at- The initial 
conditions then read 

G^/^{at, a,; k) = G,/^(0, 0; k) = ^ jj , (4.7) 



GxlA^t-. 0; k) = Gx7</'(0' k) 



1 / yM0\2 



1 - \ («.^^' 



(4i 



We introduce separate initial "temperatures" T-^^^ for the x and the modes, so 
that several scenarios can be explored. The initial masses M^^^(O) are obtained 
by solving the gap equations (13.101) at t = with = 0o- These equations are 
divergent and we use a simple local mass counterterm to take care of this. A 
much more sophisticated renormalisation scheme exists [321 [33l [MJ [351 IM], but 
this more straightforward approach has been seen to perform well [23]. For the 
coupling constant we take = \ throughout. The (renormalized) masses at 
zero temperature are m^/M^\ = 0.1 and m^/Mpi = 1. Finally, the spatial and 
temporal lattice spacings are taken as a^Mpi = 1 and at/ as = 0.01. The spatial 
lattice discretisation is fairly coarse, especially at early times when 0(t) is still 
large. The time discretisation is well under control throughout. 

We compare a range of approximations described above and in the Appendix: 

1) free: the mean field equation is solved without back-reaction, the x propa- 
gator is solved in this background (Sec. [2]). 

2) perturbative: solution of the perturbatively expanded mean-field equation, 
the X propagator is reconstructed in the perturbative coupling expansion 
(Appendix |A^. 

3) Hartree: solution of the 2PI equations keeping diagram a) (Appendix I A. 31) . 
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4) no pippi: solution to the 2PI equations keeping diagrams a) and c). 

5) full: solution of the 2PI equations keeping diagrams a), b) and c). 

Approximation 4) is used since it mimics the case considered in Ref. |37j . 
where the mean field equation includes the (perturbatively expanded) Hartree 
term, but with a dressed x propagator. Including the basketball diagram c) 
amounts to dressing the x ^"^^ V propagators with the sunset self-energy; this 
affects the mean field only via the equal-time x propagator. Comparing 4) and 5) 
also gives the possibility to assess whether a given effect originates from diagram 
b) or c). 

We now explore different scenarios and vary the initial "temperatures" Ty^j^^. 
Explicitly, we take vacuum initial conditions = 0), partly thermal condi- 

tions {T-^/ip 7^ 0, T^/x = 0), and thermal conditions {Ty./^ ^ 0). 

4.1 Vacuum initial conditions 

We initialise the system with vacuum propagators in both the x ^^nd the >^ sectors, 
= = 0. This is the relevant initial condition if we expect inflation to 
have strongly diluted all matter in the Universe, and all memory of the initial 
conditions have been lost. In order to heat the Universe, we require reheating 
(after inflation) or dissipation effects (during inflation) to create particles. 

In Fig. [3] (top) the mean field evolution in time is presented, starting in all 
cases from 0o/^pi = 10 with vacuum propagators. Note that the plot only 
shows the evolution when 0(t)/Mpi < 4 (or m^t > 7); for larger values of 
the slow-roll approximation is very accurate. The evolution in various approxi- 
mations is shown: slow-roll, free, and perturbative, and the 2PI approximations 
corresponding to the diagrams of Fig. [2) diagram a) (Hartree approximation), di- 
agrams a) and c) (no pippi), and diagrams a), b) and c) (full). The slow-roll and 
free result agree initially, but start to deviate somewhat when m^t ~ 7. Slow- 
roll ends around ^ ~ a/2Mpi, or m^t ~ 10. When including back-reaction, the 
three approximations a), a)+c), and a)+b)+c) are practically indistinguishable, 
but they differ from the free case. This is due to the time-dependent effective 
mass = + (yf^(x^(t, x)), which increases in time and causes the mean 
field to roll down faster. The time-dependent part of the mean field mass is 
shown in Fig. [3] (bottom) in the form of the subtracted equal-time (x^) corre- 
lator, (x^(t,x)) — (x^(0,x)). Here, all approximations behave in the same way, 
and approach the result predicted by the adiabatic approximation, indicated by 
the horizontal dotted line. The odd one out is the perturbative approximation, 
which can only be trusted for times shorter than, say m<^t ~ 2. This also explains 
why the mean field evolution in the perturbative approximation is much closer 
to the free one; the back-reaction on the mean field is too small. 

We find that the evolution is so slow that little departure from adiabaticity is 
observed and very few particles are created. It is therefore perhaps not surprising 
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slow-roll 
free 

perturbative 
a) 

a)+c) 
a)+b)+c) 





Figure 3: Evolution of the mean field (top) and the equal-time (x^) correlator 
(bottom) in the case of vacuum initial conditions, = = 0. The various 
approximations are slow-roll, free, perturbative, diagram a) (Hartree approxima- 
tion), diagrams a) and c) (no pippi), and diagrams a), b) and c) (full) of Fig. [2j 
The evolution in the latter three cannot be distinguished. In the bottom figure 
the dotted line indicates the asymptotic result in the adiabatic approximation. 
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that no dissipative effects are seen, and that the Hartree approximation catches 
all the main features. 

4.2 Thermal initial conditions for the x field 

We now consider the scenario where the x propagator is initially at finite tem- 
perature, T^/Mpi = 4, but the ip propagator remains in vacuum, = 0. This 
mimics the case when particle creation has been going on for a while in the x 
propagator, but has not yet trickled down into the ip propagator. In this context, 
the if propagator modes can be thought of as light particle degrees of freedom into 
which the heavy x excitations, with ~ g4'{t), may decay. This is motivated 
by the scenario where additional species are coupled to the x sector [SSI EH]. If 
dissipative effects are important, these should show up via the nonlocal diagrams 
in the propagator and the mean field equation of motion. 
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In an expanding background, the particles would be redshifted away during 
inflation, unless sourced by particle creation from the inflaton rolling. By ignoring 
the expansion in the propagator equation of motion, we keep the initial and any 
subsequently produced particles, favouring dissipative effects. We note, that as 
the X particles are very heavy initially, the initial particle numbers are very small, 
nl ~ l/[exp (MJT^) - 1] ~ 0.09. 

In Fig. Hlwe show the mean fleld evolution (top) and (x^) correlator (bottom) 
for this case. The results are very similar to the ones discussed above, although 
the back-reaction is stronger due to the x heatbath, causing the mean fleld to 
roll down slightly faster. 

We flnd therefore that the particles present in the x sector influence the 
evolution of the x propagator and the mean fleld very little|§ 

4.3 Thermal initial conditions for the (p field 

In the next scenario we initialise the if propagator in thermal equilibrium, with 
T^/Mpi = 4, and keep the x sector in vacuum, = 0. This mimics the scenario 
of Ref. |38], where integrating out the ip degrees of freedom is expected to yield 
a flnite width for the x particles. In this case particle numbers are of order 

^ l/[exp {M^p/T^p) — 1] ~ 3.5, so a signiflcant heatbath is present. As long as 
the X mass is larger than the mass, one may except the heatbath to provide 
effective damping and a thermal mass. If the coupling is strong enough, the x 
sector may thermalise as well. At the end of inflation, when M^{t) < m^, the ip 
particles in the heatbath may potentially decay directly into x excitations. 

The results for these initial conditions are shown in Fig. [5l The most striking 
effect is the change to the (x^) correlator. The presence of the ip particles af- 
fects the evolution in those approximations that are sensitive to scattering in the 
heatbath, i.e. the approximations that contain the nonlocal diagrams b) (pippi) 
and c) (basketball). The Hartree approximation does not incorporate decay or 
scattering and cannot be applied to this part of the evolution. The fact that both 
approximations a)-|-c) and a)+b)+c) show the growing (x^) correlator indicates 
that it is mostly due to the presence of the basketball diagram. We note that the 
effect is largest after slow roll ends, but is visible during the latter part of the 
slow-roll regime as well. The mean fleld evolution is very close to the one with 
vacuum initial conditions, indicating that the (p heatbath does not signiflcantly 
backreact on the mean fleld evolution using the x modes as intermediaries. It 
should be noted that the perturbative approximation does not know about the 
ip propagator, and so is identical to the case with vacuum initial conditions. 

^In Ref. [39] the x propagator was initialised in vacuum. We have repeated the simulations 
in that setup with thermal initial conditions [T-^/Mpi = 4) and found no effect of the decay 
channel as well. 
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Figure 5: As in Fig. [3l for partly thermal initial conditions, = 0, T^/Mpi = 4. 



4.4 Thermal initial conditions for both fields 

As the final case, both fields are initialized in thermal equilibrium, with T^/Mp\ = 
T^/Mpi = 4. This should further enhance the effects of the nonlocal diagrams 
and mimics the case when interactions are strong enough that the whole system 
is in equilibrium. 

The results are shown in Fig. Eland are very similar to the case just discussed 
above. The heatbath of ip particles affects the x modes mostly after the end of 
slow roll, where energy transfer from the if to x sector is effective. On the other 
hand, the heatbath of x particles makes the effective mean field mass larger than 
in vacuum, which causes to roll down slightly faster. The nonlocal diagrams b) 
and c) become essential towards the end of infiation. 

For completeness, we note that full equilibration and thermalization takes 
place on much longer time scales. This has been studied in detail using the loop 
expansion we used here in Ref . ^01 [23l HI] . Again the inclusion of diagrams b) 
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Figure 6: As in Fig. [3l for thermal initial conditions, T^/Mp\ = T^/Mpi = 4. 



and c) are essential for this. 

5 Conclusions 

The main goal of this work was to apply techniques of out-of-equilibrium quan- 
tum field theory to the problem of interacting quantum fields during and after 
a slow-roll regime. Although this topic has a long history, a number of assump- 
tions, approximations and ansatze are usually invoked along the way. Here we 
have employed a systematic formulation based on the loop expansion of the 2PI 
effective action. We considered a model where the slow-rolling field if is coupled 
to a second scalar field x via a (f'^x^ interaction. In order to separate the question 
of dissipation and particle creation during slow-roll from the dilution in rapidly 
expanding spacetimes relevant for inflation, we have treated the quantum modes 
in Minkowski space but preserved the Hubble friction term in the mean field 
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equation of motion. In this setup the mean field rolls slowly at the early stages, 
after which a smooth transition to the reheating phase is made. It also allowed us 
to bypass technical issues regarding the numerical solution and renormalisation 
in expanding space-times. 

We have studied several scenarios, corresponding to vacuum and thermal ini- 
tial conditions in the and/or x sector. The feature dominating the dynamics 
arises from the interplay between the mean field and the equal-time (x^) correla- 
tor. During the slow-roll regime, the (x^) correlator increases significantly which 
results in an increasing effective mean field mass. This effect is not important 
during the slow-roll stage, but eventually the increasing mean field mass drives 
the evolution away from slow-roll, such that this regime ends earlier than in the 
absence of interactions. We found that this is the main source of back-reaction on 
the mean field. It is already included at the level of the Hartree approximation 
and is well reproduced by the adiabatic approximation. The presence of the heat- 
baths affects the later stages of the evolution in various ways. A heatbath of x 
particles mainly speeds up the end of slow roll, since it results in a larger effective 
mean field mass. A heatbath of ip particles is important for the evolution of the x 
modes, due to the possibility of creating x particles from particles. This effect 
is important mostly at the end of and after slow roll. In this case it is essential to 
include nonlocal diagrams beyond the mean field/Hartree approximation, which 
fails to capture this effect. In all scenarios the perturbative approximation was 
found to not perform very well and to be at best qualitatively valid for very early 
times, even in the slow-roll regime. 

Particle creation and thermal back-reaction effects are crucial in warm infla- 
tion. We believe that the scenario would beneflt from being recast in the 2PI 
formalism, since in that case no quasi-particle or equilibrium assumptions would 
have to be imposed on the propagators. Also, the 2PI expansion systematically 
includes all the scattering and dissipation effects relied upon in warm inflation. 
For the parameter values used in this work, we did not reach a warm inflation 
regime. 

2PI resummed equations of motion have been successfully applied to reheating 
dynamics ^21 [T3] . There is still the issue of the numerical implementation in 
expanding backgrounds, addressed in Refs. [l3l HH HS]. Clearly, a complete 
description of inflationary dynamics requires this to be resolved. 
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A Further approximations 



In this Appendix, we present a set of approaches for approximating the local 
correlator G^{x,x) = {x^{x)). These are all well-known in the literature and are 
used in the main body of the paper for comparison. 



A.l WKB ansatz, adiabaticity and particle creation 

We consider the evolution of the mode functions, determined by 

/k(t)+c.2(t)/j,(t)=o, (A.l^ 

with 



(A.2) 



As a first approximation it is instructive to introduce a WKB ansatz of the form 

m 

/k(^) = ... exp 



— z 



(A.3) 



It is assumed that fik(^) > for all t. Inserting this into Eq. (lA.ip one finds that 



In the adiabatic limit time derivatives are small and fik(t) ~ co'k(t), such that 



"k - '^k - ^ -— + ■■■ 



2 cuk 4 cu? 



Therefore the adiabatic approximation is valid provided 



<1, 



CJk 



< 1. 



(A.5) 



(A.6) 



The (adiabatic) particle number is then given by 



nk(t) 



l + 2< 



n 



I + O [ul/ui) , (A.7) 



where is the initial x particle number present in mode k. This is the usual 
adiabatic result: particle creation is controlled by the rate of change Wk/o;^ ^6] . 

However, the absence (or suppression) of particle creation does not mean that 
equal-time correlators are constant. For the equal-time correlator 



(x(^,x)x(t 



, y)) = G^(t, t; X - y) = ^ e^-^'^^^^G^d, t; k). 



(A.8) 
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we find from tlie expressions above tfiat 

t; k) = (1 + 2<) \Mt)\' ^ (1 + 2<) ^-i^. (A.9) 

Tlierefore we find that this equal-time correlator grows in time, with a time 
dependence directly determined by the time dependence of the mean field, since 
^]<.{t) ~ (pit) for small k. Similarly, the (unrenormalized) energy density in the x 
modes, 

EAt) = / ^ (2-k(t) + 1) - / ^ (2< + 1) , (A.IO) 
decreases in time as (j){t) rolls down. 



A. 2 Exact solution 

In certain mean field backgrounds, the mode equation (lA.ip can be solved in 
closed form. In Ref. [39], we considered exponential time dependence ~ 
g-7*/2 ^Qj~ §rn?{t) ~ e~"'*) and found the analytical solution in terms of Bessel 
functionslfl We now perform a similar calculation for a slow-rolling mean field in 
a (j)"^ potential, where (j){t) is given by Eq. (12. 5p . 
It is useful to define 



and write 



/k(t) = exp ( ig(Pot + ) h^(t). (A.12) 



We find that h]f_{t) satisfies 

hk{t) + 2i6m{t)hkit) + {ul - m^) /ik(t) = 0, (A. 13) 



where = ^/k^ + rn'^ [such that uj^(t) = uj^ + 6m'^{t)]. Changing variables to 

- = (A.14) 

yields Hermite's equation, 

h'{^{x) - 2xh'Ax) + 2ukhi,{x) = 0, (A. 15) 

where 



^This solution is also relevant for a slow-rolling mean field in a </> potential. 
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The solution reads 



/ik(x) = Akif.,(x) + BkiFi ( "Y^'^' 



(A.17) 



where Hy{x) is a Hermite function and i-Fi(q!, /3, 7) is Kummer's confluent hyper- 
geometric function |17] . The coefficients ylk and i?k are determined by the initial 
conditions, /k(0) = l/-\/2^k and /k(0) = — «i^k/k(0), where f2k = ^J^l^ + 5m? {Q). 
We find that 



where 



a 



(11) 



(12) 



a 



(21) 



(22) 



A. 



VSHk^k 

-1 



(22) , -o U 
«k + "ik«k 



(12) 



a 



(21) 



zr^kftk 



(11) 



1 2 
2 '2'^° ' ' 



(A.18) 
(A.19) 

(A.20) 
(A.21) 

(A.22) 



^6-00 li^i Y' ^' ^0) ~ 2i/kV^KXo iFi - y , ^, xl^ , (A.23) 



with Xq = x(0) = 5m(0)/v^fi: and D]^ 



(11) (22) 
"k "k 



(12) (21) 
«k «k • 



A. 3 Hartree approximation 

In the Hartree approximation, the connected three-point function is neglected 
and only one- and two-point functions are preserved (Gaussian approximation) 
|48] . In that case the back- reaction of the x modes on the mean field manifests 
itself time-dependent mass and we find 

4>{t) + 3H{t)^{t) + Mj(t)0(t) = 0, (A.24) 

with 

Ml{t)=ml + g'{x\x)) = ml + g' I G^(t,t;k). (A.25) 

Jk 

As is well known |18] , the Hartree approximation as discussed here is the general- 
ization to time- dependent systems of the standard approach based on the one-loop 
effective potential in equilibrium. In thermal equilibrium, after integrating out 
the X field, one finds that 




(A.26) 
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X X 

~ + + 

Figure 7: Graphical representation of the dynamics of the x propagator (thick 
hne) in the Hartree approximation. The thin hne is the lowest order x propagator, 
the thick dashed line the propagator, and the crosses the mean field background. 
The second line shows the perturbative expansion discussed in Section [A. 4[ 



where a;„ = 27inT {n G Z) are the Matsubara frequencies. One finds that 



(A.27) 



where in equilibrium {x^{x)) is time-independent and, after performing the Mat- 
subara sum, given by 



{X\x)) 



1 + 2nB(^k) 



Here nsioj) = l/(exp(u;/T) — 1) is the Bose distribution. 

Out of equilibrium, the evolution of the quantum fluctuations can be deter- 
mined self-consistently in the Gaussian approximation in terms of the ip and x 
propagators as 



[d^ + k' + Mj(t)] t'; k) = 6c{t - t'), 

[d^ + e + M'^it)] G^{t, t'- k) = 5e{t - t'), 



with the effective masses 



Ml{t)=ml + g' 



(P\t)+ jjG^it.t- k) 



(A.28) 
(A.29) 



(A.30) 
(A.31) 



The Hartree approximation is self-consistent and easily solvable numerically. 



A. 4 Perturbative approximation 

We can also attempt to include the effect of the time-dependent mean field on 
the evolution of in a perturbative manner. This approach is indicated in 
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Fig. [71 while in the Hartree approximation all (jP insertions are summed nonper- 
turbatively, in the perturbative setup only a single 0^ insertion is preservedjj To 
obtain this, we write 

G^{x,y) ^ {Tcx{x)x{y)) = G^^\x,y) + g'G^\x,y) + 0{g'), (A.32) 

and solve Eq. (]A.29p formally order by order in g^. In order to avoid large 
corrections due to the initial large expectation value of the mean field, we write 

^\x) = 4>\0) + [4>\x) - 02(0)] =^l + A<p\x), (A.33) 

and treat g'^A(f)'^{x) as the perturbation. The expansion is therefore expected to 
be valid for early times only. We also ignore fluctuations of ip, since G^{x, x) <^ </)q 
after renormalization. In terms of the free inverse propagator, 

G^\x, y) = i[n, + ml + g^(t>l\ 5c{x - y), (A.34) 

we find the series of equations 

j^(fzG^\x,z)Gf{z,y)=5c{x-y), (A.35) 

d''zG^\x,z)G^^\z,y) = -iA<p\x)Gf{x,y), (A.36) 

etc. The formal solution of the second equation can be written as (see Fig. [7]) 

G^^\x,y) = -I j^d^zGf{x,z)M\z)Gf{z,y). (A.37) 

To first order in this expansion the mean field equation of motion then reads [37] 

[dl + mm + Mlit)] <P{t) = ig' j d'zGfix, z)A<p\z)Gf{z, x)0(x), 

^ (A.38) 

with 

Ml{t)=ml + g'Gf{x,x). (A.39) 

To make the time dependence explicit, we write the propagator along the Schwinger- 
Keldysh contour in terms of the statistical and spectral components (13.81) and find 

[df + ?,H{t)dt + Mj(t)] m = f dt'K{t,t')A^\t')<P{t), (A.40) 

Jo 

with the kernel 

K{t,t') = 2g'J d^xFf (t,t';x)pj')(t,t';x). (A.41) 
''This approximation is motivated by Refs. [49l [50l [37] . 
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In order to solve Eq. (lA.40p numerically, we need to choose a propagator for the 
memory kernel (]A.4ip . We take the simplest perturbative ansatz 



Ff (t,t';k) = !^]L±lZ^cos[cDk(t-t')], (A.42) 

pf(M-;k) = "°l"-<'-'''l , (A.43) 
^ Wk 

depending on the initial energies u)k = (k^ + + g^cfl)^ . One may think of more 
elaborate ansatze for these propagators and include e.g. more thermal effects as 
well as the effect of scattering and decay processes [161 EZ] • However, we prefer 
to discuss this in a systematic approach and employ the 2PI effective action to 
consider extensions beyond the mean field approximation. 

Finally, one can consider further approximations to these equations, using a 
derivative expansion [l9l [50], [37]. However, we will not do so here. 
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